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^-^, Abstract 

The Sachs equations governing the evolution of the optical matrix of geodetic WANDs (Weyl 
aligned null directions) are explicitly solved in n-dimensions in several cases which are of interest in 
potential applications. This is then used to study Einstein spacetimes of type III and N in the higher 
dimensional Newman-Penrose formalism, considering both Kundt and expanding (possibly twisting) 
solutions. In particular, the general dependence of the metric and of the Weyl tensor on an affine 
parameter r is obtained in a closed form. This allows us to characterize the peeling behaviour of the 
^ . Weyl "physical" components for large values of r, and thus to discuss, e.g., how the presence of twist 

affects polarization modes, and qualitative differences between four and higher dimensions. Further, 
5— ( , the r-dependence of certain non-zero scalar curvature invariants of expanding spacetimes is used 

^Wi to demonstrate that curvature singularities may generically be present. As an illustration, several 

explicit type N/III spacetimes that solve Einstein's vacuum equations (with a possible cosmological 
constant) in higher dimensions are finally presented. 

> 

(^ ! 1 Introduction 

m 

^Nj ' It was recognized long ago that several important features of gravitational radiation in General Rel- 

ativity can be conveniently described in a covariant manner by studying asymptotic properties of 
spacetimes [1-4]. More generally, the development of asymptotic techniques has proven fundamen- 

^^ ' tal in understanding general properties of the theory, since the behaviour of the gravitational field 

near (spacelike or null) infinity encodes essential information about physical quantities such as mass, 
angular momentum and flux of radiation (at least in asymptotically flat spacetimes). From a tech- 
nical viewpoint, the Newman-Penrose formalism [3] turns out to be extremely useful in analyzing 

^^ , the fall-off properties of gravitational fields at infinity. In a nutshell, (after making certain initial 

}_^ ' technical assumptions) from the Ricci identities one first extracts the r-dependence of some of the 

Cw ' Ricci rotation coefficients which are needed in the analysis (throughout the paper r will denote an 

affine parameter along null geodesies). Subsequently, specific Bianchi identities are integrated that 
determine the behaviour of the Riemann (Weyl) tensor, which geometrically characterizes properties 
of the gravitational field. From this one demonstrates, e.g., the characteristic peeling-off properties 
of radiative spacetimes [1-4] . Then one can further proceed to integrate the remaining Ricci/Bianchi 
identities and thus find asymptotic solutions (possibly with specific extra assumptions, see, e.g., [5], 
and [6] for a review and further references). 

It was noticed in [1] that four dimensional algebraically special spacetimes, while leading to sig- 
nificant mathematical simplification, still asymptotically retain the essential features of (outgoing) 
radiation fields generated by more realistic sources. In particular, in that case the r-dependence of 
the Weyl tensor can be determined in closed form (and not only asymptotically). This is also an 
important first step towards the exact integration of the full Newman-Penrose equations, aimed at 
determining the explicit metric functions. One can indeed get quite far in the case of algebraically 
special spacetimes, at least in vacuum [7-13], and several important exact solutions fall within this 
large class (see [6,14] for further references). 

In recent years the interest in gravity in higher dimensions has grown considerably, mainly mo- 
tivated by modern unified theories, AdS/CFT and recent brane world scenarios. Notions such as 
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the total energy of an isolated system and energy flux are thus fundamental also in higher dimen- 
sional theories [15-17]. The study of radiation in higher dimensions may ultimately enable one to 
distinguish different models, since properties of gravitational waves depend on the model under con- 
sideration (in particular, on spacetime dimensions [18]). It is thus now of interest to explore ideas and 
techniques similar to those mentioned above for the study of asymptotic properties of algebraically 
special spacetimes with n > 4 dimensions. The necessary formalism has been provided in recent 
works, where an n > 4 generalization of the Petrov classification [19] and of the Newman-Penrose 
equations [20-22] have been presented.^ Such a programme was thus started in [24], where we studied 
vacuum spacetimes admitting a non-degenerate (as defined in section 2) geodetic multiple WAND, 
and satisfying a further condition necessary for asymptotic flatness. Thanks to the results of [20], 
the assumed non-degeneracy implies that the only possible algebraic type of that family is II (or D) , 
but not III and N. We observed that such spacetimes do not peel-off and do not contain gravitational 
radiation, as opposed to the case of four dimensions [1]. 

The analysis of [24] can now be extended in various directions by modifying some of the assump- 
tions made there. It is the purpose of the present paper to focus on empty spacetimes of type III 
and N, with a possible cosmological constant. For these, the (unique) multiple WAND is necessarily 
geodetic and degenerate [20]. We can already remark at this point that this implies that such space- 
times can not be asymptotically flat (even in a "weak", local sense, in contrast to the n — A case).^ 
Our analysis can be however still of interest for spacetimes with, e.g., Kaluza-Klein asymptotics. It 
is also worth observing that, as opposed to [24], we will not need here any extra assumptions on the 
asymptotics of the Weyl tensor - its full r-dependence will be fixed by the Bianchi identities.^ 

The paper is organized as follows. In section 2 the Sachs equations for a congruence of geodetic 
WANDs £ are studied for spacetimes that satisfy the (rather weak) condition _Roo = Rab^''^'' ~ 0, 
and an explicit solution is given when the principal directions of shear and twist are aligned. This 
includes several cases of interest and, in particular, Einstein spacetimes of type N and III, which are 
then studied in the rest of the paper. Namely, in sections 3 and 4 we determine the r-dependence 
of the Weyl tensor components of such spacetimes in a parallelly transported frame. We can thus 
discuss, e.g., their peeling behaviour near infinity, curvature invariants and possible singularities, 
frame freedom and rotation of frames induced by the presence of twist. Differences with respect to 
the four-dimensional case are also pointed out. In the final appendix we construct several examples 
of Einstein spacetimes of type N and III. These are explicitly given in five dimensions, but they can 
also be easily extended to higher dimensions if desired. 

Notation Following [19-22], we use a frame consisting of two null vectors m(o) = i (aligned 
with the multiple WAND) and m^i^ = n, and n — 2 orthonormal spacelike vectors mj^), where 
i,j,... = 2,...,n — l. In terms of these, the metric reads 

Qab = 2l(anb) + 5rjm^a^m[^\ (1) 

where, hereafter, a, b — 0, 1, . . . , n — 1. 

Derivatives along the frame vectors i, n and m^i) are denoted by D, A and 5i, respectively. We 
choose the frame such that it is parallelly transported along I. The optical matrix L of .^ has matrix 
elements 

Lij — (.a;bm.1i)m(j-j, (2) 

with (anti-)symmetric parts 

Sij — L(ij^, Aij = L[ijj. (3) 

The optical scalars expansion, 9, shear, a, and twist, to, are defined hy 6 = Lii/{n — 2), a^ — 
{Sij — 05ij){Sij — OSij), and oj^ — Aij Aij. Other Ricci rotation coefficients used in this paper are 
defined by (see [20,22] for the full set of coefficients) 

Lii — la-.bn^rn^i), Ln = ia-.bmfi^n , Ln = £a;bn'^n , Mjk — —Mik = 'm^i)a;b'm'lj)m(^ky 

(4) 
The Weyl tensor of spacetimes of type III and N has only negative boost weight frame components, 
for which we use the compact symbols 

"^i ~ Cioii, ^ijfc = —Cikij, 'iij ~ -^C\i\j. (5) 



^Very recently, an extension of the GHP formalism to higher dimensions has also been developed [23]. Although we 
will not need the GHP formalism here, the results of [23] are useful also in the NP context, since some redundancy of the 
original Bianchi equations of [20] has been removed. Beware of the fact that some normalizations used in [23] slightly difEcr 
from those of the present paper. 

2 Of. footnote 1 of [24]. 

^The extra assumption needed in [24] concerned the asymptotic behaviour of the Weyl components Cijki (in the notation 
of [19], sec also the following), which vanish identically in the case of type III/N spacetimes. 



From the symmetries and the tracelessness of the Weyl tensor one has the identities [20] 

*i = 2*jjj, *{ijfe} = 0, *ijfe = -*jifc, *ij = *3i, *ii = 0, (6) 

where *I'{ijfc} = ^^fe + ^kij + ^ jki, which will be employed throughout the paper. A subscript or 
superscript will denote quantities which are independent of r (e.g., s^j), ao, ^Ijk^ etc.). 

2 Geodetic WAND: solving the Sachs equations 

We are interested in studying asymptotic properties of spacetimes along the congruence generated 
by a geodetic multiple WAND £. As a first step, it is thus natural to fix the r-dependence of the 
matrix Lij , which determines the optical properties of i. This can be done by integrating the Sachs 
equations (a subset of the Ricci identities), which for a geodetic, affinely parametrized WAND (not 
necessarily multiple) read [22] 

OLij ^ — Lik-^kj: (7) 

where we assumed that the Ricci tensor satisfies Rqq — (which obviously holds for Einstein spaces, 
defined by Rat ~ QabR/n). When Lij = we have a trivial solution of (7) corresponding to Kundt 
spacetimes. In the rest of this section we will study only the non-trivial case Lij 7^ 0. 

When £ is non-degenerate, i.e. L is invertiblc, such a matrix differential equation can be easily 
solved in terms of L~^ [24], and by taking the inverse matrix one then finds L. When the number of 
dimensions is kept arbitrary, this is done more conveniently by expanding L in a power series in 1/r. 
This was indeed the starting point in the analysis of [24]. By contrast, here we will not assume L 
to be invertible, but we will solve (7) in a closed form under some other assumptions. These will be 
such to include the form of L compatible, in particular, with type III/N spacetimes [20]. However, 
in this section we will be general enough so that the presented results will apply in a wider context. 

2.1 Explicit solution when the principal directions of shear and twist 
are aligned 

If one forgets for a moment about the request that the frame be parallelly propagated, one can always 
choose the basis vectors mj^) such that the symmetric or the antisymmetric parts of the matrix L 
take their canonical form (that is, Sij is diagonal, or Aij is block-diagonal with 2-dimensional anti- 
symmetric blocks). It is natural to refer to such preferred basis vectors mj^) as the principal shear 
directions and principal twist directions, respectively (there may be some degeneracy, in general, i.e. 
there need not be a unique basis of principal directions). Here we consider the special case when 
there exists a basis of vectors that are principal directions of shear and twist simultaneously, which is 
relevant to important applications (and includes, in particular, the case when L is a normal matrix, 
i.e., [S, A] = 0/ This means that L admits a canonical form given by a direct sum of 2- blocks of the 
form 

C^f 1^'f^O, (8) 



^ -^23 S(3) , 

where, for definiteness, the frame indices refer to the first block (next blocks will be characterized by 
pairs of indices (4, 5), (6, 7), ... ). If the spacetime dimension n is odd, there will be also an extra 
one-dimensional block. 

2.1.1 The canonical frame can be parallelly transported 

We want to solve eqs. (7), which hold when the basis vectors m(i) are parallelly transported along I. 
Since we are now assuming that the mj^) coincide with the common principal shear and twist direc- 
tions, for consistency it is necessary to prove that such vectors can indeed be parallelly transported. 
Let us thus take such a canonical frame at a special value of r, say r — 0. Then we have, by 
construction, that I/jr=o is block-diagonal, with blocks of the form 

^l'-°=f ^.'^'fO- (9) 

V -^23 S(3) J 

One can now define a frame in a neighborhood of r = by parallel transporting the frame defined 
at r = 0. By Taylor-expanding L = L\r=o + r{DL)\r=o + |r^(D^L)|r=o + . . ., the evolution of L 
(assumed to be analytic) will be determined by (7) and its r-derivatives 

£>™L = (-l)'"m!L'"+\ (10) 

evaluated at r = 0. Since L"^\r^o has clearly the same block-diagonal structure (for any m), it follows 
that {D"^ L)\r=Q has the same block- form as well. Consequently L has such a block-diagonal form for 
any r, with blocks given by (8), that is what we wanted to prove. 



2.1.2 Explicit r-dependence 

Since now each 2-diniensional block obeys a decoupled equation, we have effectively reduced eq. (7) 
to the standard Sachs equations of n = 4 gravity. For the first block, we thus have the general 
solution [1,3,25] 



sip) +r[s^2)S(3) + 



{A'2S?] 



l + K^(2)+^(3))+^'K)4) + (^°: 



23) 



(P = 2,3), (11) 



A" 

A23 = p r '-^, ^, (12) 



and similarly for the next blocks. In odd dimensions, the last block will be 1-dimensional and given 

by S(„_i) = S('„_i)/(1 + rs^n-i))- 
The optical scalars are 



More special subcases are discussed below. 

2.1.3 Non-twisting case 

The case oj = is of course included in the above "aligned" case since Aij — 0. Then, Lij — Sij is 
diagonal with eigenvalues (using (11) with A^j = 0) 

^ + rs(i) 

as obtained in [26]. 

Note that in this case £ must be expanding {6 7^ 0) [22] or one is simply left with the Kundt class 
(see, e.g., [27-29] and references therein). 

2.1.4 Non-shearing case 

When a — one has Sij — 6Sij, which is clearly aligned with any Aij. Then S(i) — 6 for any 
i — 2, . . . , n — 2, so that by (11) all non-zero Aij take the same value (up to a sign). All the 
information about Lij is thus contained in the two optical scalars 



(^o + r{9i + ^^u^i) 



LOO 



1 + 2reo +T^{el + :^Lon 1 + IrOo + r^ 91 + -^^^ul 



(15) 



This special case was already discussed in [22]. Similarly as in the non-twisting case, we necessarily 
have here 9 =^ [22] (one can also easily see that detL 7^ and £ is thus non-degenerate), unless we 
consider Kundt spacetimes. If one has, in addition, also cj = 0, one is led to the class of Robinson- 
Trautman spacetimes [30]. In fact, the case oj 7^ is possible only if n is even [22], so that in 
odd-dimensions non-shearing spacetimes thus reduce to either the Kundt or the Robinson- Trautman 
class. 

2.2 Case of a matrix L satisfying the "optical constraint" 

A special case of a matrix L admitting aligned principal shear and twist directions arises when L 
satisfies the so called "optical constraint" [31], i.e. (dropping the matrix indices) 

[S,A]^0, A^^S'^-TS, (16) 

where J-" can be fixed by taking the trace (note, in particular, that L is thus a normal matrix). 
Considering matrices L satisfying this special property is motivated by the fact that this includes 
several important classes of vacuum solutions, namely spacetimes of type III and N [20], Kerr-Schild 
metrics [31] and general asymptotically flat spacetimes with a multiple WAND [24]. In particular, 
in four dimensions the optical constraint is equivalent to the shearfree property (except when ^ = 
and, simultaneously, rank(5') = 1). 



One can easily see that under the above assumptions the matrix Lij takes the block diagonal form 
(cf. also [31]) 



^(1) 












£-(p) 



\ 




(17) 



V 

The first p blocks are 2x2 and the last block £ is a (n — 2 — 2p) x (n — 2 — 2p)-dimensional diagonal 
matrix. They are given by (after appropriately rescaling r, see [31] for details) 



^M - 



S{2m) 



S(2m) 

r 



^2m,2m+1 
S(2m) 



(m = i,---,p), 



(a 






£ = -diag(l,...,l, 0,. 
r 



2fi,2fi+l 



«{2m) 



^•^ + K2.))" 



(18) 

(19) 



(m — 2p) (n — 2 — m) 



with Q<2p<m<n—2. The integer m denotes the rank of Lij, so that L^j is non-degenerate when 
m = n — 2. 

For certain purposes it turns out to be very convenient to define the complex combination 



P(M) 



S(2m) 



iA 



2/j,2/j+l — 



"(Zm) 



This satisfies the compact Sachs equation 



Dp(f,) = -p{p), 



(20) 



(21) 



which generalizes the standard Sachs equation of the four-dimensional theory [1,3,14,25] (except that 
our p(^) differs by a sign from the standard NP symbol). 

In the rest of the paper we will study spacetimes for which there is at most one non-zero block (for 
^ = 1, say). In such a situation we can drop block-indices and introduce the more compact notation 
s = 8(2) ~ L22 = L33, A = A23 = L23 = —L32, p = P(i) = s -I- lA. If s = then also ^ = [22], 
i.e., one has the trivial solution L = corresponding to Kundt spacetimes. On the other hand, for 
non-zero expansion s 7^ one finds 



s + iA- 



- tao 



(22) 



which determines the expansion scalar 6 — 2s/ (n — 2) and the twist oj = \/2A of £. 

Let us emphasize that all the above results hold for a geodetic WAND, not necessarily a multiple 
one, under the only assumption Roo = on the energy- momentum content of the spacetime. In 
the following we will however consider more special situations and in particular restrict to geodetic 
multiple WANDs of Einstein spacetimes. 



3 Type N spacetimes 



In the previous section we have determined the r-dependence of the optical matrix Lij for a wide class 
of spacetimes. This matrix plays a special role since it enters the Bianchi identities that need to be 
integrated in order to find the r-dependence of the Weyl tensor. The full set of Bianchi identities has 
been presented in [20] (cf. also [23]). In the following we will use those containing D-derivatives, i.e., 
derivatives in the direction of the geodetic multiple WAND i. We start our analysis by considering 
type N Einstein spacetimes. The case of type III spacetimes is similar but technically more involved, 
and will be dealt with in the next section. 

According to [20], one has to consider the two possible cases Lij — and Lij 7^ (the results 
of [20], relying on the Bianchi identities, were obtained in the case A = 0, but hold also for Einstein 
spacetimes since, for these, Rabcd;e ~ Cabcd;e ~ the same applies in the next section in the type III 
case) . 



3.1 Kundt spacetimes 

The condition Lij = defines Kundt spacetimes (i.e., £ is non-expanding, shearfree and twistfree). 
In this case the Bianchi eq. (23) of [20] (derived from (B.4, [20])) reduces to 

£i*y = 0, (23) 

so that 

*.. = *?i, (24) 

does not depend on r. The amplitude of the gravitational field is unchanged as one moves along light 
rays (the multiple WAND), which is a "plane wave" -like behavior. By performing a space rotation, 
one can thus always align one's frame to the "polarization axes", i.e., to the eigendirections of ^^j. 

We note in addition that Einstein spacetimes of type N belonging to the Kundt class are VSI 
(if A = 0) or CSI (if A 7^ 0) spacetimes (i.e., spacetimes for which all scalar invariants constructed 
from the Riemann tensor and its covariant derivatives are either vanishing or constant) [21,32]^ and 
no physically useful information can thus be extracted from their invariants. Finally, it may be also 
worth recalling the result of [31] that, in vacuum (with A = 0), Kundt solutions of type N coincide 
with non-expanding Kerr-Schild spacetimes. 

3.2 Expanding spacetimes 

When Lij does not vanish identically the results of [20] implies that L is normal and in particular 
of the special form obeying the "optical constraint" discussed in section 2.2. In an adapted frame 
(which we take to be parallelly transported, cf. section 2.1.1) it has the only non-zero components 
s = S(2) — S(3) = (n — 2)6/2, and A23 = —A32 — A, as given in (22). 

Further, it follows from [20] that in such a frame "^ij has the only non-zero components'' 

*33 = -*22, *23=*32. (25) 

Bianchi equation (23) in [20] can thus be compactly written as 

D(*22+i*23) = -p(*22+i*23), (26) 

with the simple solution 

*22 + i*23 = (*22 + **23)P- (27) 

Under a rotation in the m'^'-m''^' plane such that the frame is still parallelly transported along £ 
we have ^22 + i'^23 — >■ e'"''(^22 + **I'23), hence the ^23 (or >I'22) component can be set to zero if and 
only if ao = 0, i.e. if the twist vanishes. In that case one has simply »l/22 ~ ^22/1' (or ^23 ~ ^23/^) 
and such frame is thus aligned with the eigenframe of ^ij (cf. also [26]) or rotated by 45°. In the 
twisting case this is not possible and the eigendirections of ^ij will spin with respect to the parallelly 
propagated frame, i.e., the effect of twist is to "mix" the two polarizations as one proceeds along the 
rays of the gravitational field (indeed, in four dimensions ^22 and ^23 correspond to the real and 
imaginary part of ^4, respectively, i.e. the well-known "-I-" and "x" polarization modes). 

Eq. (27) can be employed, for instance, to characterize the asymptotic behavior of the Weyl 
tensor (whose leading term clearly behaves as 1/r, as in four dimensions). In addition, knowing 
the r-dependence of the Weyl tensor can also be useful for studying possible spacetime singularities. 
Namely, one can analyze the simplest non-trivial curvature invariant admitted by type N spacetimes 
in four [34] and higher [21] dimensions, i.e., 

r — /^aibia2b2;ciC2 y^ ^ei die2d2 1/1/2 /^ { OQ\ 

1M = '~J Oaidia2d2;ciC2'-^ ^ei6ie2 ^2 ;/l/2 • 1^°/ 

It can be shown® that for type N Einstein spacetimes 7jv is proportional (via a numerical constant) 
to 

r r^,T, ^2 , /,T, ^2^2/ 2 , .2x4 (*22) + (*23) ,„„. 

/iVOc[(*22) +(*23)J {S +A) = . 2 , „2N6 " (29) 



*The "VSI part" of this statement has been proven in [21] in the case A = 0. When A 7^ (i.e., for the "CSI part") the 
proof goes essentially unchanged, since the Bianchi identities and most of the Ricci identities which one needs are unaffected 
by the cosmological term. Even though for A 7^ Nij is not a balanced scalar, Lemma 4 in [21] still holds and the Weyl 
(but not the Riemann) tensor and its derivatives are also balanced and thus all its invariants of all orders as well as mixed 
invariants with the Ricci tensor vanish. The only difference from the vacuum case is that one can now construct non-zero 
constant invariants using contractions of the Ricci tensor R^b = dabR/i^ (see also a remark at the end of Section 1 of [33] 
in the four-dimensional case). A similar comment will apply also later in the case of type III Einstein spacetimes of the 
Kundt class and will not be repeated there. 

""We observe that in the canonical form given in [20] one has >I'23 = 0. However, the frame used in that paper was not, 
in general, parallelly transported. We will comment more on this point shortly. 

®This was done in [21] for vacuum spacetimes with A = 0, but it can be easily extended to the case A ^ 0. 



If ao (which generically is a function of coordinates other than r) vanishes at some spacetime points, 
then there will be a curvature singularity at r = = ao- In particular, this always occur in the 
non- twisting case [26]. Further ("r-independent") singularities may arise from a possible singular 
behaviour of the function (^['22)^ + ('I'm)^ (see appendix A.l for specific examples). Note also that 
/jv — > for r — > 00, i.e. far away along the multiple WAND. 

4 Type III spacetimes 

Similarly as for type N, also in type III vacuum spacetimes the matrix Lij has, in an adapted parallelly 
propagated frame, the only non-zero components s = S(2) = S{3) — {n ~ 2)6/2 and A23 = —A32 ~ A 
[20],^ given in (22) (the fact that an adapted frame can be taken to be parallelly transported follows 
from section 2.1.1, as in the type N case). 

In order to fully determine the r-dependence of the Weyl tensor we shall integrate the Bianchi 
equations (B.l), (B.4), (B.6) and (B.9) of [20], which contain D-derivatives of negative boost weight 
Weyl components. These equations (given later in an appropriate context) contain, in addition to 

i 

Lij, also the Ricci rotation coefficients Ln, Ln and Mjk, and the derivative operator St. We thus 
first need to fix the r-dependence of these quantities. Together with Ln these will also determine 
the behaviour of the metric coefficients. 

4.1 Ricci identities 

The assumption that £ is the multiple WAND of type III spacetimes and that the frame is parallelly 
transported greatly simplifies the general form of the Ricci identities given in [22] . Recalling that the 
r-dependence of L is given in (22), we can use the results given in Appendix D of [31].* It will also 
be convenient to divide the space indices i,j, k into two groups, corresponding to the non- vanishing 
and vanishing block of Lij, respectively, i.e. 

(30) 



(31) 
(32) 

(33) 

r + hi, (34) 

where, hereafter, for brevity lowercase symbols In, In, rriji^ and In denote quantities independent of 

i i j 

r. Because of the index symmetries of M J fc [20], we require m^fe -t-TTj^fe = for any i, j, fc = 2, . . . , n — 1. 

4.2 Commutators and r-dependence of the metric 

In order to determine the radial dependence of the derivative operators Si and A let us take the afBne 
parameter r as one of the coordinates and x any set of (n — 1) scalar functions (which need not be 
further specified for our purposes) such that (r, x ) is a well-behaved coordinate system. Then the 
directional derivatives (when acting on scalars) take the form 

D = dr, A = Udr + X^dA, 5^ = cj,9, + itdA, (35) 

where Oa = d/dx^. 

The r-dependence of these can be determined using the following commutators [21] 

AD - DA = LiiD + LaSi, (36) 

6iD - DS, = LuD + Lj,5j. (37) 

Applying (37) on r we get Dui — —Lu — LjiUjj, which, using (22) and (31), leads to 

uj2 + ii^3 ^ {i^2+i^3)p^ {h2 + ilri)rp, ujn, = —hwr + uj^. (38) 



p,q,o = 2,3, v,w,z — 


4,5, ...,n- 1. 


The relevant coefficients thus read 


Li2+iLr.i ^ {Ii2+ili3)p, Li^ ^ h^. 


L21 + iLii — {hi + il3i)p, Lnii — lull. 


i i i i i i 

Mj2 + iMj3 = {'mj2 + imj3 )p, Mjw = nij^ , 


ill = Rc[(;i2 + ill3){l21 - il3l)p\ + 


R , , 

n{n — 1) 



^More precisely, this has been proven in full generality for all n > 4 non-twisting solutions of type III and for all ?i = 5 
solutions of type III (in Appendix C of [20]). For twisting type III solutions with n > 5 it was assumed that a "generality" 
condition on the components ^ijk holds and that ^I"; ^ (Section 4 of [20]). We do not consider possible "exceptional" 
cases in this paper. (Note also that the possibility rank(S) = 1 was not discussed explicitly in section 4 of [20]. However, 
it can be shown from equations (69), (71), (74)-(77) and (80) therein that such case indeed cannot occur.) 

®Up to a small difference: Appendix D of [31] studied Ricci-flat spacetimes of type II (or more special), while in this 
section we consider Einstein spacetimes of type III. 



Similarly, acting with (37) on x gives D^^ = —Lji^j , so that 

6 + «C3 = (C2 + <3 )P, (.w =£.n, ■ (39) 

Applying (36) on x leads to _DX — — La^i , from which (using (39)) 

X^ = Re[(/2i - ihi){£.°2 + iCaV] - ln,iCr + X^° . (40) 

Applying (36) on r gives DL'' — — Ln — Lnuji. Using also (32), (34) and (38) we obtain 

R 



t/ = Re 



[a;2 +ia;3 -r(Zi2 ~ ihz)](l2i -ikn)p\ 



2n{n - 1) 



^Ik; ^k 



r'-(«ii+/„ia;^)r + t/°. (41) 



Note that the above expressions for the Ricci coefficients and the derivative operators have been 
derived under the assumption p 7^ (i.e., Lij 7^ 0). When this does not hold, there is no need to 
distinguish between two types of indices as in (30) and the corresponding expressions can be obtained 
from the above results simply by dropping all quantities containing indices 2 or 3. 

Let us observe at this point that, since in the above coordinates 

r = 5^, n" = U5^ + X'^S'X, mj',) = uj,5^ + ^^5%, (42) 

the r-dependence of all components of the frame vectors is now known. This automatically also fixes 
the radial dependence of the metric (1). For the contravariant components we explicitly have 

rr nrr , rA -^^A , /-A AB ^A^B r AO\ 

together with (41), (38), (40) and (39). The covariant components can be found by imposing the 
orthonormality relations among the frame vectors, which gives 

grr = 0, X^grA = 1, it QvA = 0, 

2(7 + X^'X^ffAs = 0, ui,+X'^£,fgAB=Q, itif9AB = 0. (44) 

The explicit form of the covariant coefficients can be worked out more conveniently after further 
information about the line-element is specified and, possibly, other adapted coordinates are defined, 
and we will not discuss this any further here (see, e.g., [26] for the non-twisting case). 

Wo now proceed with determining the r-dependence of the Weyl tensor. As in the type N case, 
let us discuss the two possible cases Lij = and Lij ^ separately. 

4.3 Kundt spacetimes 

When Lij = the Bianchi equations (B.6) (or, equivalently, (B.9)) and (B.4) of [20] take the form 
D*«jfc = 0, (45) 

k 

2D*ij = 5j^i + *i(Lij - Ljx) + 2^jkiLki + *fc Mij ■ (46) 

Direct integration of the first of these gives 

*..fc = *?,*, *, = *? = 2*0,.,.. (47) 

We now discuss eq. (46). Since p = for Kundt spacetimes there is not need to introduce two 
types of indices and for the Ricci coefficients we have simply 

i i 

L\i = l\i, Lii = In, Mjk ~ ™3fc , (48) 

while the derivative operator St reads 

5,^{^lur + Lj'^)dr + if°dA. (49) 

Using (47), (48), and (49) the integration of (46) gives 

*,, = i [if'i'lA + *"(«!, - hi) + 2*;'fc.«fci + -^Im), ] r + *°,.. (50) 

Furthermore, since ^a = = '^[ij] the integration constants appearing above must satisfy 

*">=0, C«^"*?,A + *?(?!, -2/a+r«;,)=0, (51) 

*fo] = 0, C*°i,A + *° («|i|.] - ^.]i) + -^'hkhi + -^Iml^] = 0. (52) 

The above equations (47) and (50) thus fully describe the r-dependence of the Weyl tensor for 
type III Einstein spacetimes of the Kundt class, in agreement with [26] (where an adapted frame such 
that hi — hi ~ was used). Note that, in contrast to the type N case, now Weyl components of 
boost weight —2 do in general depend on r. This is the typical peeling-off of Kundt spacetimes (here 
restricted to type III) and is well-known also in four dimensions [1]. As discussed for the type N, also 
Einstein metrics of type III that belong to the Kundt family fall in the VSI or CSI class [21, 32]. 



4.4 Expanding spacetimes 

Bianchi equations (B.l), (B.9), (B.6) and (B.4) of [20] read 

L»*, = -2*fcLfe„ (53) 

D-^jki = "S/iAjk + -^kuLij - <S>jiiLik, (54) 

2D*yfe = -*iLjfc + *jl/ifc - 2<i'ijiLik, (55) 

2D*ij = -2^ikLkj + 5j*j + "ifiiLij - Lji) + 2*jfciifci + *fc My . (56) 

We now study the above differential equations for various index combinations (recall (30)). 

4.4.1 Components of boost weight 1 

In terms of the two index sets (30), eq. (53) can be conveniently rewritten as D{^2+i'i^3) = ~2p{^2 + 
1*3) and D*u, = 0, so that 

*2 + i*3 = (*2 + i*3)p', *» = *»• (57) 

Let us consider (54) and (55) for the components "i/wvp- These reduce, respectively, to D\E'iu„p = 
and -D(\E'mi,2 +i'^mv3) = —p{'^wv2 +i'i^mv3), whose only common solution is clearly (since p 7^ here) 

*».p = 0. (58) 

Before integrating the remaining Bianchi equations containing the operator D, we take advantage 
of the fact that for type III spacetimes some other Bianchi identities become purely algebraic. These 
are eqs. (B.7), (B.ll) and (B.16) of [20], where a detailed analysis can be found. In particular, from 
(B.7) and (B.16) one can derive eq. (58, [20]), which reads 

e{n - 2)*,,fc + 45[,|,*,fc|,j - 2Ssk'i'^,s + 2S'[,|fc*|,] = 0. (59) 

Since 9 j^ 0, using this with {i,j,k} — {v,w,z}, {2, 3, w}, {p,w,v}, {w,p,q} we get, respectively 
(recall also "i/iijk} ~ 0) 

VE'i.u.z = 0, '$23w = *»32 - *»23 = 0, '^p^v =0, f^u = 0. (60) 

In addition, from ^i = 2^ijj the last two equations give 

*233 = -*2, *322 = -*3, *»33 = -*u,22. (61) 

Since the r-dependence of \l/23p is thus now determined by \&p, the only remaining equation to be 
solved is eq. (55) for *uj22 {~ — 'I'tu33) and '^ui23- This can be written as 

D(*»22 + 4*^23) = -p(*»22 + i*»23), (62) 

with solution 

*™22 + i*»23 = P (*^22 + 4*^23) • (63) 

This fixes the r-dependence of all boost weight —1 Wcyl components for type III Einstein space- 
times. One can check that all boost weight +1 and Bianchi equations given in [23] (thus, in 
particular, cqs. (53)-(55) and the above mentioned algebraic equations) are now satisfied. 

4.4.2 Constraints on the Ricci rotation coefficients 

Before we proceed with fixing the r-dependence of the boost weight —2 Weyl components, it turns 
out that suitable Ricci identities will lead to considerable simplifications useful in the following cal- 
culations. Let us consider Ricci equation (Ilk, [22]), which for a geodetic null congruence £ in type 
III Einstein spacetimes reduces to 

^lj\^i\k] = Li[j\Li\k] + LiiLyk] + LiiM[jk] + Li[j\Mi\k]- (64) 

Considering the various equations obtained for i,j = p,q and k = w and using (22), (31)-(33), 
(35), (38) and (39), one finds 

n^22 = "133 =i^ui, "i23 =— m32 =— /i«,ao— 4u, ao,A- (65) 

Next, for i, j, k = o,p,q one gets 

«(C2 +^£,3 )ao,A = —{uJ2 +i'^3) +iao[—{li2 +ili3) + 2{l2i +il3i)], (66) 



and for i = w^ j = q and k = z 

2 3 

m^z =0 = m„^ . (67) 

Finally, for i = w, j = 2 and fc = 3 we find m^a = Imiao so that, by (65), 

m23 =-ms2 =lwiao, ^w'^ao,A = -aoQwi + hm)- (68) 

Other index combinations do not contain any further information. We also note that all the 
remaining Ricci identities [22] contain also Ricci rotation coefficients that do not appear explicitly in 
the Weyl tensor components and thus we do not consider those here. 

4.4.3 Components of boost weight 2 

As the last step, the r-dependence of boost-weight —2 components of the Weyl tensor can now be 
determined from eq. (56). Using the above results for ^i and 'i'ijk, it is convenient to study various 
cases with different values of the indices. 

For i,j — w,z, recalling (33) and (67), eq. (56) becomes simply 

I>*^„. = 0, (69) 

so that 

*... = *L, *[».]= 0, (70) 

where the latter condition follows from '^wz ~ 'i'zw 

Next, for i = w and j = 2,3 eq. (56) can be written as 

21>(*i„2 + i*»3) = -2(1/„2 + i*i„3)p + *2(Mi„2 + *M„3) + *3(Af»2 + iM^a)- (71) 

Using (57), (33), (65) and (68), this leads to 

*»2+i*»3=p(*^2 + i*™3)+^°p', P» =i(a;° +iaoi„i)(*5 + i*^), (72) 

with *{'„2] + «*[»3l = 0. 

We also observe that for i — w, j — 2,3 the antisymmetric part of (56) gives (since ^Pjij] = 0) 

S^i'i'2 + i*3) = V [{hw - Li + iml^ )(*2 + i*3) + 2(*L + J*»3) + 2(*L2 + **°23)a2i - ihi)] 

- 2p\loI + ia„/„i)(*':^ + i¥i,), (73) 

from which (with (35), (38), (39) and (68)) 

- it^i-^l + i*3),A = (*2 + i*3)(3;i» - ln,l + irn-i^ ) + 2{hl - i/3l)(*»22 + **;i23) 

+ 2(<2 + **°3)- (74) 

Finally, we have to consider (56) in the case i,j — 2,3. The corresponding equations can be 
compactly rearranged as 

2D(*22 + *33) - (<52 - i<53)(*2 + J^s) = (1'2 + «*3) [(il2 - ^^13) - 2(L21 - iLai) 

+ (M33 + JM32)] - 2P(1'22 + *33), (75) 

2D(*22 - *33 + 2^*23) - (^2 + i5:i){'^2 + «*3) = (*2 + 1*3) [(^12 + 1^13) + (-M33 + iM32)] 

- 4L„l(*»22 + i*»23) - 2p(1'22 - *33 + 2i*23)- (76) 

Using (35), (38), (39) and (57) we find the needed transverse derivatives, i.e., 

(52 + i<53)(*2 + i*3) = 2p^(*'2' + i¥l)[i{(i" + i^3")ao,A - (4' + ^4') + r{li2 + ilia)] 

+ p'(6''° + Je3''°)(*2+i*3),A, (77) 

{52 ~ i<53)(*2 + i*3) = 2p^p(*J; + i*^)[i(C2'"' - i^f)ao,A - (a;2" - i4') + r{li2 - iha)] 

+ p'p(6^"-iC3''")(*2+i*3),A. (78) 

After substituting (77) into (76), by direct integration we find 

*22 - *33 + 2i*23 = p-''^" + P^B" + pC" + prV\ (79) 
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where (using also (66)) 

^0 = (*o + j*0) [{uj° + iJi) - iaoihi + ihi)] , 

S" = -i [(*2 + i*3)[3(/l2 + i/13) + ^(m32 + im's )] + i^f' + ie3'"')(*2 + i*3).A] , 

C" = *'2'2-*«3+2J*^3, (80) 

Note that here C'' is the only new integration "constant" . 
Similarly, by substituting (78) into (75), one has 

*22 + *33=Pp:^°+pg°, (81) 

where 
J^ = -\ [(*2 + i*3)[3(/i2 - ilvi) + (m33 + Jm32 ) - 2{hi - ikn)] + {^f ~ i?;f )(*2 + i*").^] , 

e" = *!^2+*33, (82) 

and Q'^ is the new integration constant. 

However, since "^u = 0, we have *22 + ^33 ~ —'^mw = — *S,ii,- This must now be compatible 
with (81), from which we get 

*L=0, *22 + *33=0, F° = Q, (83) 

so that 

*22+*33=0 = *„„. (84) 

AAA Summary and discussion 

The above results can now be conveniently summarized as follows. First, a number of Wcyl compo- 
nents vanish identically, namely 

*» = 0, *ij„ = 0, ■^n,z2 = = *„^3, (85) 

*22 + *33 = 0. (86) 

The r-dependence of the non-zero components is given by 

*2 + i*3 = 2(*233 + i*322) = p'(*2 + **3), (87) 

*™22 + 1*^23 = --^wii + i'^w32 = p(*S,22 + «*»23), (88) 

*... = *L (*[».] = = *'^„), (89) 

^^2 + i-^^Z = p{^l2 + **L) + p'Vl, (*|'„2] = = *|'„3]) (90) 

2(*22 + J*23) = P^^" + p'b" + pC° + prV'', (91) 

where the various integration "constants" satisfy (72), (80), and (83) with (82), along with the 
constraints (66), (68) and (74). 

As r — > oo, one can observe a different fall-off behavior for different components and therefore 
a peeling- like behavior. There exist components of boost weight —1 both with l/r^ (eq. 87) and 
1/r (eq. 88) leading terms. The slower fall-off described by the latter equation can be qualitatively 
understood as due to the fact that there is no expansion along the "w-directions" . As for boost weight 
—2, in general there are components that are asymptotically constant in r (eqs. (89) and (91)) and 
components that fall off as 1/r (eq. (90)). Again, the asymptotically "constant" terms are due to the 
non-expanding extra-dimensions. 

There are several special subcases that may be worth mentioning. First, for the special subtype 
Ill(a), which is invariantly defined by the condition Vl'i = [19], one obtains the simplifying conditions 
po = ^0 ^ 23" = and, by (74), VE'S;,2 + 1*^,3 = -(^21 - «?3i)(*^22 + **»23) (further simplification 
can be achieved by using a residual frame freedom, see below). Next, also in the non-twisting case 
(p = 1/r) the above expressions (87)-(91) become much simpler and were given already in [26] (in 
particular, thanks to (66) one gets ^ — 0, so that the p^ term of (91) disappears). Finally, one 
can compare the above results with the well-known asymptotic behaviour in four dimensions (cf., 
e.g., [6, 12, 13]). In that case, in our notation, indices v, w, z do not exist (since i, j, fc = 2, 3 only) and 
eqs. (87) and (91) encode all the information, corresponding, respectively, to the complex Newman- 
Penrose scalars ^tg ~ 1/r^ and ^4 ~ 1/r (note that I?" = in four dimensions). There are no 
non-expanding extra-dimensions and thus no terms with a slower fall-off. 

In the general case, for certain applications the asymptotic behavior of the Weyl components (87)- 
(91) can in fact be visualized more clearly by taking a series expansion. Using p — X]m=i(*'*o)™~^r~™, 
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p2 = Em=i™-(««o)'""V-('"+^) andp^ = 5E^=i'^("T'+l)(*ao)'""V-('"+2), up to the leading and 
sub-leading terms one finds 



T^Zr 'V'O / \ T* 



*»22+J*»23= (^^^ ;3— J + ^ (^^— + ,.2 J + g'(^' )' (93) 

Wu,2 + «Wu,3 = 1 ^ ^^-^ 1 + * 1 \- 



2|.2 J \ r 2r^ 



+ 0{r-''), (94) 



r,0 , *22 +"0^tul*'i23A , ,f , ,t,0 , *23 " «0^«;1*„22 \ , ^nr„-2\ 



*22 + i*23 = ( -/»1*L2 + ^^ ' ™'-^M + i ( -U*L3 + '■' — ---^M + 0(r-^), 

(95) 

while still ^l/u,^ = Vl'SLz- This clearly demonstrates, in particular, how the presence of twist mixes 
up various polarization modes. Simpler expressions can be obtained by performing specific frame 
transformations, as briefly discussed below. 

Similarly as for the type N, the r-dependence of the Weyl tensor can also be used to discuss the 
possible presence of curvature singularities. The simplest non-trivial curvature invariant for expanding 
type III Einstein spacetimes is [21,35] 

r /^aibia2b2;ei ^ /^dici d2C2 ;e2 ^ tQf<\ 

J^III — ^ ^aicia2C2;ei v^ ^di6id2&2 ;e2 • \^^) 

It can be shown that [21] 

IlII ex (S^ + A^f[9{^^<il^f + 27(*,\E'i)(*«.22*u,22 + *»23*»23) + 28(*,„22*»22 + *»23*u,23)^] 

' [i^'ir + i^^^f]' , ^^ mr + (^ir](^l22^i22+^l23^l23) , ,,,^,0 ^0 ^^0 ^0 .2' 
y r-5— 2T2 1" ^' 2~, 2 <' 20(,Wu,22Wti,22 + W«,23Vm23J 

(r +aQ) r +aQ 



1 



(r^+ag) 



2-14 • 



(97) 



As in the type N case, there may be curvature singularities localized at points where r^ + Oq = 0, 
which may or may not exist, in general (but they always do in the non-twisting case [26]). Addi- 
tional singularities may also arise from a possible singular behaviour of \E'2, ^3, Vl/Jl,22 and ^$^23 (see 
appendix A. 2 for specific examples). 

4.4.5 Prame freedom 

The results above have been obtained using a generic parallelly transported frame and therefore hold 
in any such frame. For certain purposes it may be desirable to simplify some expressions by using 
the freedom to perform null rotations 

£^ £, n^ n + z,m}-^ - -ZkZk£, m! -^ m<'^ - zd, (98) 

where Dzi = in order for the new frame to be still parallelly transported. 

For instance, once can set to zero the Ricci rotation coefficients I/12 + iLiz or L21 + iI/31 (by 
taking Z2 + izs = —{I12 + ilia) or Z2 + izs = — (/21 -I- i^3i), respectively), or one may want to set 
to zero certain Weyl components of boost weight —2 (for type III spacetimes components of boost 
weight —1 are invariant under (98)). Namely, one can transform away the term B'^ in (91) by taking 
2(2:2 -I- i23)(*2 + 4*3) = 13° (note that if *2 + **3 = then B° is automatically zero and there is no 
need to perform any null rotations). Alternatively, if ^2,22 + *'J'il)23 7^ one can set \E''^2 + *^2j3 in 
(90) to zero by taking (22 — i23)(^m22 + ^"^^23) ~ '^m2 + i-'^ms- Additionally, there are null rotations 
in the ui-directions. If *I'S,22 +*^ti)23 7^ one can take 4:Zw{'iw22 +'i'^w23) ~ — C° so as to have, in the 
new frame, C" = in (91) (note that this is not possible in four dimensions since ^^22 + **I'il;23 ~ 
in that case). Or one can set P^ = in (90) by taking 2t„(*2 + «*§) = 2P^ (if ^% -|- i*§ = then 
V'ij = already in the original frame). 

Furthermore, spatial rotations can also be used to simplify the form of certain Weyl components. 
For example, a spatial rotation in the m(2)-m(3) plane adds an arbitrary (r-independent) phase to all 
the above non-zero components (except for ^P^z, which is unchanged) and can thus be used to set to 
zero the imaginary part of the corresponding integration constants (in particular, in the non-twisting 
case p is real and one can thus align one's frame to the "polarization" of such components). Next, 
one can use rotations in the planes defined by the "non-expanding" directions Tn(„) to, e.g., align the 
frame to the direction defined by *I'uj22 or ^u)23, etc., or to diagonalize '^^z- The most convenient 
way how to use the frame freedom may depend on the specific spacetime under consideration and its 
possible symmetries. 
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5 Concluding remarks 

After presenting some general results about the Sachs equations (section 2), we studied specific 
features of Einstein spacetimes of type N and III in arbitrary higher dimensions. This is a natural 
extension of previous studies such as [20,26] and partly complements, in different respects, other works 
either by the present authors or by others, e.g., [24,29]. In particular, by explicitly determining the 
r-dependence of the Weyl tensor we were able to discuss several physical properties of the general 
families of solutions of type N/III, either with or without expansion, and to compare these with 
their well-known four-dimensional counterparts. The results of this paper also represent a first step 
towards the exact integration of the full Newman-Penrose or GHP equations for such spacetimes, 
which will be studied elsewhere. In the following appendix the discussed results are illustrated by 
presenting some explicit solutions that, to our knowledge, have not been given before. 
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A Some explicit expanding spacetimes 

In the main text we have studied properties of general Einstein spacetimes of type N and III in higher 
dimensions, for an arbitrary value of the cosmological constant. While Kundt solutions are similar to 
their four-dimensional counterparts and several explicit examples are already known [27-29,31,32,36], 
not many type N/III spacetimes with Lij ^ have been found. To our knowledge, in fact, the only 
such examples have been obtained (for A = 0)^ as a direct product of a four-dimensional type N/III 
Ricci-fiat spacetime with an Euclidean Ricci-flat space (see, e.g., [26]). The lack of less trivial examples 
is partly due to the fact that, contrary to the four-dimensional case, they are necessarily shearing [20] 
and therefore they do not show up, e.g., in the Robinson- Trautman family [30]. 

In this appendix we present a few examples of such solutions (both non-twisting and twisting) 
which are not direct products. They are in fact warped products and solve the vacuum Einstein 
equations Rat ~ 2A.gab/{n — 2) with a possible cosmological constant (which can take an arbitrary 
value, at least in the non- twisting examples). As a reader familiar with four-dimensional exact 
solutions may easily note, these spacetimes have been constructed by appropriate "warping" of four 
dimensional type N/III solutions (see, e.g., sections 13.3.3, 28.1, 28.4, 29.1-29.4 of [14] and section 19.2 
of [38]; some of the original papers are also quoted below in the appropriate context). 

In fact, all the considered metrics can be written in the form (cf [39]) 

ds' = j^dz^ + /(^)da^ (A.l) 



where 



2 . „ , , . , 2A 



/(.) = -A.^ + 2d. + 6, A=^— ^^^— ^, (A.2) 

b and d are constant parameters,^'' and Aa^ is a Lorentzian Einstein spacetime of dimension n ~ 1 
with Ricci scalar 

R„^{n-l){n-2){\h + d'). (A.3) 

This metric will be specified in the following and will characterize the properties of the full spacetime 
ds^. We observe that the latter is a direct product only in the special case of a constant f{z), i.e. 
\ — Q = d (with & > 0). In order to have a Lorentzian signature for ds^, we require f{z) > 0, which 
may restrict possible parameter values and (possibly) the range of z. Namely, since f{z) has real 
roots if and only if R„ > 0, when Re < we require A < {R„ = admits also A = 0, but this 
case simply corresponds to a direct product), while for Re > any sign of A (including A = 0) is 
admitted, at least for suitable values of z. 

A.l Einstein spacetimes of type N 

A. 1.1 Non- twisting case 

Using the above general ansatz (A.l), one can obtain five-dimensional type N Einstein spacetimes with 
an arbitrary value of the cosmological constant A by taking da^ to be the general four-dimensional 



^An Einstein space which is the direct product of non-Ricci-fiat Einstein spaces also contains Weyl components of boost 
weight [37] and thus can not be of type N/lIl. 

'^'^In fact, A is the only physically meaningful free parameter contained in f(z) (as one can always redefine z ^ az + l3) 
but for convenience we will generally keep also b and d unspecified. 
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expanding and non-twisting type N Einstein metric with a possibly non-zero (four-dimensional) Ricci 
scalar Ra = 12{\b + d^). This was given in [40] (see also [38] and references therein, in particular for 
a transformation to the standard Robinson- Trautman coordinates) and reads 

do-^ = -2ipdudr + 2r^(da;^ + dj/^) - 2r(2r/i + ex)dudx - 2r(2r/2 + ey)dudy + 2(^B + A)du^, (A.4) 

with 

A = \e\x^ + y") + e(/i:r + hy)r + (/f + jl)r\ (A.5) 

V = l + i6(a;2+y^), (A.7) 

where e = ±1 or and the functions /i = fi{x,y) and /2 — f2{x,y) are subject to 

d.fi = a^/2, dyfi = -a./2. (A.8) 

In the case \b + d^ — the metric da^ is Ricci-flat and the spacetime ds^ can be lifted to any 
higher dimensions by simply replacing da^ — > da^ + X]a('^-^<^)^- One can obtain a higher dimensional 
solution also in the case \b + d^ 7^ 0, however in a bit more complicated way (a simple direct product 
will not work as it will introduce Weyl components of boost weight zero, cf. footnote 9). For simplicity, 
in the following analysis we will restrict to the n = 5 case. 

The geodetic multiple WAND is given by 

£ = dr, (A.9) 

with r being an afhne parameter along the corresponding null geodesies. We can then choose a 
parallelly transported frame in the form 



1 « 

On 



^f{z 



f\b + d^ _\\ 2_ e/2 + rd^] 1 / ex\ 

[ f(z) 2)' ^f{z) J* ^f{z)V' + 2r J '''' 

1(2) = j==dx, m(3) = j===dy, m(^4) = —j==rdr + \J f{z)dz- (A.ll) 



y27(i) ^' -^ r^w^)"' -^ ^m 

Then the only non- vanishing Ricci rotation coefficients relevant to the discussion in the main text 
are (note that L\\ 7^ but we do not need it here) 

T T ^ T ^iJ -^'^i^ + W) r r Xz~d 

M2 — Las ~ -, -L21 +tJ^3i = , , , J ^14, = —J^4i = F==- (A. 12) 

r r^/2fiz)^ ^/f{z) 

The non-vanishing independent components of the Weyl tensor are (after using (A.8)) 

*.+.*. ^-%f)^i. (A.13) 

The curvature invariant In given in (28) is therefore proportional to 



[{dif2f + {dif2) 



,212 



Similarly as in the four-dimensional case [34], in five dimensions Im diverges whenever any of the 
following conditions hold: i) r = 0; ii) V = (i.e., for e — —1 and x + y = 2); iii) the quantity 
(p'xf'i) + (pyf'i) diverges. In five dimensions an additional curvature singularity is also located at 
the roots of f(z) — 0, which are present iff Ra > 0. 

A. 1.2 Twisting case 

Here we present a five-dimensional twisting Einstein spacetime of type N with a negative cosmological 
constant A. This is constructed by taking da^ in (A.l) to be the four dimensional type N twisting 
solution of Leroy [41] (but in different coordinates, cf. also [14,42,43]) with a negative Ricci scalar 
R„ = UiXb + d"^) = -4s2, i.e.. 



'^°" = -2^ 



-(r^ + l)(da;^ + dy^) + -{dx + y^du) [Qydr + j/^(l - r^)d« + (13 - r'^)dx + 12rd2y] 



(A.15) 
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The coordinate r is an affine parameter along the geodetic multiple WAND £ = dr- We choose a 
parallelly propagated frame 

n = — du +W2dr dx +widy — (Xz — d)rdz, (A. 16) 

"^(2) = ^^'(yff 1) ^ i^9u + Vra. - y\dx ~ y'dy) , (A.17) 

(a„ + V9. - 2/^9:. + j/Vaj,) , (A.18) 



^ SyS/S 1 1^ ^ ^ 2a ,3a ,,.3 



m(4) = ^^=4ra. + y7(iy9,, (A.19) 



//( 



z 



where 



3/(2:)(r2 + l) 2 6/(2)(r2 + l) 



The non- vanishing components of the optical matrix are 



1-22 + «i23 = ^33 - 4^32 = :• (A.21) 

r — I 



The remaining relevant non-zero Ricci rotation coefficients are 



1/21 - Ji31 = 2l M32 + «M33 = -xhjr^' Ll4=-L41= 7==, A.22 

r-iy3/(2) V/(2) 

2 2 3 3 _j ;^2 _ d 
M42 + iM43 = -i(M42 + «M43) = : , (A.23) 



The Weyl tensor components are 

*22 + i*23 = TT, r^-jj-T^ , (A.24) 

in agreement with the general result (27). Hence for the curvature invariant In (28) we have 

Due to above mentioned relation Re, = — 4s'^ < 0, the warp function f{z) has no real roots and 
therefore in this case Im is everywhere regular. In addition, the components of the Weyl tensor in the 
above frame (parallelly propagated along i) are also regular. See [42] for a discussion of the regularity 
of the four-dimensional Leroy metric (A. 15). 

To conclude, we note that Leroy's solution was rediscovered in [43] as a special case of a more 
general class of four-dimensional twisting Einstein spacetimes of type N (determined up to solving 
a system of two third-order ODEs). Similarly as above, these metrics can be used to construct 
other type N solutions in five dimensions (in particular, also Ricci-flat ones if one starts from a four 
dimensional geometry with a positive Ricci scalar). Finally, let us also observe that further five (or 
higher)-dimensional type N solutions with a negative cosmological constant can easily be constructed 
by taking da^ to be the well-known four dimensional type N Hauser solution [44] (cf. also [14]). 

A. 2 Einstein spacetimes of type III 
A. 2.1 Non-twisting case 

An explicit solution in n = 5 dimensions is given by (A.l) with 

do-^ = ^(da;^ + dy^) -t- 2dudr + ( ^x + (\b + d^)r^ ) d^t^ (A.26) 

a;3 y 2 J 

which is a four-dimensional Robinson- Trautman spacetime of type III [14]. The cosmological constant 
A can take an arbitrary value. Similarly as in A. 1.1, when \b + d^ = the above spacetime can be 
lifted to any higher dimensions by simply replacing da^ -^ du^-l- 5]]^(dza)'^, but again in the following 
we will restrict to the n — 5 case. 

A geodetic multiple WAND is given by 

e = dr, (A.27) 
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while a suitable parallelly transported frame consists of 






\b + (f _ \\ 2 33a; 
f{z) ^2)'' +32/(z) 



3x^ 



4/(z)r 



dx - {\z - d)rdz 



(A.28) 



1(2) 






„3/2 



1(3) 



-dy 



The only non-zero components of the optical matrix Lij are given by 



m,(4) = —^===rdr + \JT{z)dz 



(A.29) 



The remaining non-zero Ricci rotation coefficients relevant to the discussion in the main text are 



L\2 — L21 — M33 



3^/S 1 



For the non-zero Weyl tensor components one finds 

3x^/^ 1 



Li4 — —Lai — — 



\z — d 



V^ 



*2 = 2*233 = 



4/3/2(z)r2' 



*24 = 



?,x-"\\z-d) 1 



8/2 (^) 



(A.30) 



(A.31) 



which are special cases of the general expressions (87) and (90). 
For the curvature invariant (96) we now have 

IlII (X 



,P(z)r^ 



(A.32) 



There are thus obvious curvature singularities at r = and a; — > cx) and, in the case \b + d^ > 0, 
additional curvature singularities are present at the roots of f{z). 

Other four-dimensional Robinson- Trautman spacetimes of type III [14] can be also used instead 
of (A. 26) to obtain different five-dimensional solutions. 



A. 2. 2 Twisting case 

A five-dimensional solution with a negative cosmological constant is given by 



/(^) 



-\z' 



(A<0), 



(A.33) 



and by taking da^ to be the four-dimensional Ricci-flat type III metric [14] 



da" 



^ + 6y xw I da; H -dy + 2(\/l3 + l)yx w dxdy 

'^3 / ,^3 



— xwdu i 6ydx + (\/l3 + l)a;dy J — 4ywdrda; + 2dMdr + -xdu 
where the functions w = w{x) and S = S(r, x) are defined as 



(^/T3-5)/2 



■^ 2,62 

X, — r + X w . 



(A.34) 



(A.35) 



The constant parameter a gives rise to non-zero twist, and for a = (i.e., kj = and E = r^) the 
solution corresponding to (A.34) reduces to the previous (A. 26) (in the case b — = d). 
As above, the coordinate r is an affine parameter along the geodetic multiple WAND 



e^dr 



(A.36) 



A parallelly transported frame is now given by 



1 /, 3rwa; w . „ 



A 2 33a; 

"2'" ^ 32/(2) 



3(H-V13)a:'«; 



8S/(2) 



3ra; 

47Rs' 



3x''w 



\2) 



Hs) 



1(4) 



^3/2 



2ryxwdu + ( |s + i(l -|^ ^/l3)a;^wM dr + xrd^ + x^wdy 



\z 



X w —2yx'wdu + -(1 + V 13)r9r — xdx + rdy 



(A.37) 



TTR 



rdr + \/fiz)dz 
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The only non-zero components of the optical matrix Lij are given by 

1 

r + ix'^w 



1/22 + iL2:i = L33 - iL-i2 = ZTTTZTT.' (A. 38) 



Other relevant non-zero Ricci rotation coefficients are 

.,2 . 2 , 3Jx 1 , ^ Az 

L12 -I- ^Ll3 = L21 - i-Lai = -«(M32 + iMss) = ^== — , . o , ii4 = -L^i = 



2 .2 3 .3 Aza;^™ i , , . 

M42 + «M43 = -« M42 + «M43 = —r^= , . 3 • A.39 

y/f{z) r + ix^w 

The non-zero Weyl tensor components are given by 

3a;*/2 1 

*2+i*3 = 2(*233+i*322) = -TT7irTr^— 1-T^, (A.40) 

*22 + 1*23 = -^^-777727^^ f , . 3 ^3 > 1'24 + tl'34 = o .2 ^ ^ ^ i ■ 3 ^2 ' ^.41 

16/2(2:) (r + ja;^w)^ 8/2(2:) (r + «a::'^w)2 

These can be compared with the general expressions (87), (91) and (90). 
For the curvature invariant (96) we thus have 

x^ 

^"' °^ \6 12/ 2 — T^Te • (A.42) 

^b_3i2^j..2 j^ rj.b^2jb 

Curvature singularities are thus located at r = for a; — > (more in general, the regularity at 
r = — X depends on how this location is approached) and at 2: = (while for a; — >■ cx) we have 
Iiii — >■ 0, as long as a 7^ 0). 

Note that the above five-dimensional spacetime can be lifted to any higher dimensions by replacing 
da^^da^' + J^JdZc.f. 

It is finally worth observing that a four-dimensional twisting spacetimes of type III with a negative 
cosmological constant was given in [42] (see also [14]), while several Ricci-flat solutions are known, 
see [14] and references therein. These can all be used instead of (A. 34) to obtain other five-dimensional 
solutions. 
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